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The calculation of the 0(a s ) virtual corrections to the matrix element of the in- 
clusive decay b — » s-y is reported. These contributions drastically reduce the large 
renormalization scale dependence of the leading logarithmic calculation. Com- 
bining these results with the preliminary result for the Wilson coefficient C7(mt,) 
calculated recently by Chetyrkin, Misiak, and Miinz, we estimate the branching 
ratio to be BR(B -* X s j) = (3.25 ± 0.50) X 10" 4 . 



1 Introduction 

Rare B meson decays are particularly sensitive to physics beyond the standard 
model (SM). In order to extract the effects of possible new physics, precise 
experimental and theoretical work on these decays is required. 

The inclusive mode B — > X s ^ can be systematically analyzed with help 
of the expansion in inverse powers of the (heavy) 6-quark mass, m^. At the 
leading order in such an expansion, the inclusive decay rate is given by the 
perturbatively calculable free quark decay rate. As the power corrections start 
at the 1 / ml level only, we neglect these contributions in this talk and model the 
decay rate T(B — ► X s j) by the quark level decay width T(b — > sj), including 
perturbative QCD corrections. 

The earherSJVI computations of the branching ratio for B — > X s ^ presented 
e.g. in refs. BB'tru are in full agreement with the CLEO measurementB BR(B — > 
X s j) = (2.32 ± 0.67) x 10~ 4 . There are large uncertainties, however, in both 
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the experimental and the theoretical results. In view of the expected increase 
in the experimental precision, the calculations must be refined correspondingly 
in order to allow quantitative statements about new physics. So far, only the 
leading logarithmic QCD corrections of the form a™ log" (m^/m^,) have been 
resummcd completely A systematic improvement is obviously obtained by 
calculating all the next-to- leading terms of the form a s a™ log" (m^/m^). 

Before discussing the various steps leading to a next-to- leading result, we 
briefly review the theoretical framework in which the process b — > s~f(+g) is 
evaluated. Usually one starts form the effective Hamiltonian which is obtained 
by integrating out the i-quark and the VF-boson. Neglecting power (mb/m w ) 
suppressed operators of dimension > 6, the effective Hamiltonian reads 

AC 8 

H eff (b^sj) = ^VtbVZ^CiMOiM , (1) 

where the quantities Cj(/x) are the Wilson coefficients evaluated at the renor- 
malization scale fx and the Oi are following operators: 

02 = (c La ^b La ) (sl^Clp) , 

7 = (e/167r 2 ) s a {m b {^)R + m s {^)L) b a , 

Os = (g s /16TT 2 )s a a^(m b (^R + m s (^L)(X^/2)b p G^ u . 

As the Wilson coefficients of the penguin induced four-Fermion operators 
0^,...,Oq are very small, we do not explicitly list them here. 

From the /i-independence of the effective Hamiltonian, one can derive a 
renormalization group equation (RGE) for the Wilson coefficients Cj(/z): 

M^0»(M) =7j»0j(M) i ( 2 ) 

where the (8 x 8) matrix 7 is the anomalous dimension matrix of the operators 
Oi. Working to leading- logarithmic precision only, it turns out that it is suf- 
ficient to do the matching (at /1 = m w ) neglecting QCD corrections; to solve 
the renormalization group equation using the order a s anomalous dimension 
matrix 7^°^ ; and to calculate (perturbatively) the matrix elements of the op- 
erators Oi at the scale \i w m b , again neglecting QCD corrections. Athough 
it is clear that the renormalization scale /1 should be of the order of mj, its 
precise value is not fixed, of course. Following common practice, we vary /j, 
in the range rrib/2 < [i < 2m b . This variation leads to a large (±25%) scale 
dependence of the leading logarithmic result. Analytically, the source of the 
large scale dependence is due to a term of the form ~ a a log([i /rn%). 
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2 Steps to a next-to-leading result 



In order to get the next-to-leading logarithmic result for the branching ratio, 
one has to improve the Wilson coeffients at the scale fj, rs mi, and in addition 
one has to work out the 0(a s ) corrections to the matrix elements for b — > sj. 
The improved Wilson coefficients are obtained in two steps: First, the match- 
ing at the scale [i — m\y has to be calculated including order a s corrections!!! 
Second, the RGE step down to the scale /i w mi has to be done using the or- 
der a 2 8 anomalous dimension matrix 7^. This second step is the hardest one, 
because some entries of the anomalous dimension matrix (like 727) have to be 
extracted from 3-loop diagramsQ! The calculation of the order a s corrections 
to the matrix element b — > sj involves the bremsstrahlung process b — > sjg 
and virtual corrections to b — > sj. While the bremsstrahlung corrections (to- 
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gether with those virtual corrections whi 
singularities) have been worked out earlierl 
worked out the virtual corrections completely El. Technically, the most difficult 
part was the calculation of the order a s corrections to the contribution from 
the operator O2', the corresponding 2-loop diagrams are shown in ref. a. Using 
the Mcllin-Barnes representation of certain progagator type denominators, we 
could write the result M2 of the 2-loop diagrams in the form 

M 2 = w + £c„(^)\o g ™(<) , (3) 

with n = 1,2, 3, 4, ... and m = 0, 1, 2, 3. The coefficients cq and c nm are pure 
numbers, i.e., independent of any parameters like mb,m c , ... . Note, in partic- 
ular, that there is no naked log(m^/m^) term present in eq. (^). So the limit 
m c — + of Mi exists. 



3 Preliminary results for the branching ratio BR(B — > X s ^j) 

Summing up, in order to get the next-to-leading logarithmic result for BR(B — > 
X s ~/), one has to know both, the 0(a s ) matrix elements and the next-to-leading 
order Wilson coefficients f] at /i ps m&. Only the combination of these two in- 
gredients is independent of the renormalization scheme. It turns out that in 
the naive dimensional scheme (NDR) with MS subtraction, the correction to 
CV(wb) is smallQ. Therefore, a good approximation for BR(B — > X s ^j) is ob- 
tained by using the leading value ior Cy(mb) in the numerical evaluation of the 
matrix elements, as presented inEl. While the /1 dependence was about ±25% 
in the leading logarithmic calculation (varying /1 between m^jl and 2m;,), it 

b In fact, it is sufficient to know only Ct{i±) to next-to-leading precision. 
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Figure 1: Branching ratio for B — > X s -f as a function of rat- The upper (loweri solid curve 
is for fi = mt,/2 (/x = 2rrij,). The dotted curves show the CLEO 1 — a boundsQ. The other 
input parameters are taken at their central values. 



gets drastically reduced to ±6% when taking systematically into account the 
virtual corrections to the matrix elements. The term ~ a s log(ml/ [i 2 ), which 
caused the large scale dependence of the leading logarithmic result, is cancelled 
by the 0(a s ) virtual corrections to the matrix element. In Fig. 1 we show the 
remaining ^-dependence as a function of the top quark mass mt] all the other 
input parameters are taken at their central values. Combining the uncertain- 
ties in m t and fi (m t — (170 ± 15) GeV; mfe/2 < fx < 2mb) leads to an error 
of about 9% in the branching ratio. Besides that, there are other errors to be 
taken into account, stemming from the uncertainties in a s (mz), the semilep- 
tonic branching ratio, and the ratio m c /ro&. Taking a s (mz) = (0.117± 0.006), 
BR s i = (10.4 ± 0.4)%, and m c /m b = (0.29 ± 0.02), one obtains an extra 
error of about ±12% El. To conclude, we end up with a preliminary predic- 
tion for the branching ratio BR(B -> X sl ) = (3.25 ± 0.30 ± 0.40) x 10~ 4 , 
where the first error is due to the (fi,mt) variation and the second error 
due to the other input parameters. Adding the errors in quadrature, we get 
BR(B -> X sl ) = (3.25 ± 0.50) x 10~ 4 . 
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